Equality connecting energy dissipation with violation of fluctuation-response relation 
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In systems driven away from equilibrium, the velocity correlation function and the linear response 
function to a small perturbation force do not satisfy the fluctuation-response relation (FRR) due to 
the lack of detailed balance in contrast to equilibrium systems. In this Letter, an equality between 
an extent of the FRR violation and the rate of energy dissipation is proved for Langevin systems 
under non-equilibrium conditions. This equality enables us to calculate the rate of energy dissipation 
by quantifying the extent of the FRR violation, which can be measured experimentally. 
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Macromolecules and colloidal particles of the order of 
nanometers to sub-micrometers suspended in an aqueous 
solution provide an ideal ground to study the founda- 
tion of non-equilibrium statistical mechanics. Recent ad- 
vances in experimental techniques for the observation and 
manipulation of such small systems has generated the 
possibility of directly investigating the non-equilibrium 
nature of fluctuations in the system. In particular, such 
techniques have been designed to verify several universal 
relations such as the fluctuation theorem 0, Q, the 
Jarzynski equality 0, Qi and the Hatano-Sasa identity 
110. 

Through the investigation of non-equilibrium systems, 
it has been recognized that the quantification of the vio- 
lation of the fluctuation-response relation (FRR) || pro- 
vides new information for systems driven far from equi- 
librium On the other hand, the rate of en- 
ergy dissipation is the most fundamental quantity that 
characterizes non-equilibrium steady states. Thus, it is 
naturally expected that the FRR violation is related to 
the amount of energy dissipation. Toward this end, in 
this Letter, we present an equality between the rate of 
energy dissipation and an extent of the FRR violation for 
a class of non-equilibrium stochastic systems. 

In this Letter, for simplicity, we mainly study a system 
described by the Langevin equation: 



jx(t) = F(x(t),t)+at)+eF(t), 



(1) 



where 7 is a friction coefficient, and £(t) is zero-mean 
white Gaussian noise with variance 27T. In particu- 
lar, we consider two models of the force F(x,t). In 
the first model, which is termed Model A, F(x, t) = 
-S a (t)od x Uo(x) - S^xdxUxix), where U l {x) for i = 0, 1 
are periodic potentials with period t, and a(t) is a Pois- 
son process on {0, 1} with a constant transition rate a for 
the transitions from both to 1 and I to 0. This model 
was originally studied as a model of motor proteins, and 
is termed a flashing ratchet model \vA Il3| . In the second 
model, termed Model B, we address a time-independent 
force F(x, t) — F(x) = f — d x U(x), where / is a constant 
driving force and U(x) is a periodic potential with pe- 



riod £. See Ref. 0] for the physical significance of this 
model. The last term in Eq. with a sufficiently small e 
represents a "probe" force used to investigate the linear- 
response property of the system. The initial condition is 
given at t — t ln \t and we consider the limit i; n ;t — -> — 00 
below. 

First, we define several measurable quantities for this 
system. Assuming v s to be the steady-state velocity with 
e = 0, the response of the velocity to the probe force is 
characterized as 



(x{t)) e = v s + , 



R(t-s)F(s)ds + 0(e 2 ), (2) 



where (• • •) denotes the ensemble average in the pres- 
ence of the probe force with e. The response function 
R(t) has the causality property R(t) — for t < 0. An- 
other important quantity is the time-correlation function 
of velocity fluctuations in the absence of perturbation, 
which is defined as 



C(t) = ({x(t) - v s ] [±(0) - v s }) . 



(3) 



Next, we quantify the energy dissipation in the 
Langevin description. According to the definition in 
Ref. ^3 , the rate of energy dissipation J(t) at time t 
for each trajectory is expressed as 

J(t)dt = [ 7 ±(t) - f (i)] o dx(t), (4) 

where o denotes the Stratonovich multiplication |l6| . 
This definition of the rate of energy dissipation complies 
with both the conservation of energy and the second law 
of thermodynamics [ljj . 

Based on the above, we present a Theorem that con- 
stitutes the main claim of this Letter: 
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C(u>) - 2TR'{u) 



2^ 



(5) 



We express the Fourier transform of an arbitrary func- 
tion A(t) as A(uj) = J_ A(t) exp(iujt)dt. The prime 
denotes the real part of the function. It is widely ac- 
knowledged that in equilibrium, i.e., for cases wherein 
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FIG. 1: C(u) and 2TR'(u) as functions of w when D = 5 
and a = 10. All the quantities are dimensionless under the 
normalization asj — £ — T=l. The solid and dotted lines 
represent C(uj) and 2TR'(lj), respectively. 



a = for Model A and / = for Model B, the correlation 
function C(t) is connected to the response function R(t) 
as C(t) = TR(t) for t > 0, which is the FRR in model 
d) [EtlSl Thus, from C{t) = C{-t) and R(t) = for 
t < 0, we obtain the Fourier transform of this relation as 



C(u) = 2TR'{uj). 



(6) 



Here, we stress that Eq. © does not hold in non- 
equilibrium steady states; thus, it is observed that the 
right-hand side of Eq. J^} represents the extent of the 
FRR violation. 

It is noteworthy that the equality (JSJ) holds for variety 
of Langevin systems regardless of the magnitude of ex- 
ternal driving as well as the manner in which the system 
is driven away from equilibrium. We will provide a de- 
tailed explanation of this equality in another paper. In 
the present Letter, we demonstrate the validity of Eq. (JjjJ 
by the numerical verification for Model A and provide a 
mathematical proof for Model B. 

First, we investigate Model A. Statistical quantities 
for Model A can be calculated by analyzing the following 
Fokker-Planck e quat ion that corresponds to the Langevin 
equation [l2,0j|: 



d_ ( Pa(x,t) 
dt \Px{x,t) 



C - a 
a 



a 

C\ — a 



P (x,t) 
Pi(x,t) 



(7) 

where d = -d x [Fi(x) + e/ p (<) - Td x ]h and F^x) = 
—d x Ui(x) for i = 0, 1. In particular, the potentials were 
selected as Uo{x) = Dcos(2ttx/£) and Ui(x) = const., 
and all the quantities were converted into dimensionless 
forms by normalizing 7, £, and T to unity. The func- 
tions, C(ui) and R(ui), were calculated accurately by the 
matrix-continued-fraction method [l4|. In Fig. 2] C(u>) 




FIG. 2: Energy fluxes (J) and (Ji n ) calculated on the basis 
of Eq. and Eq. (|SJ , respectively, as functions of the tran- 
sition rate a. The solid line and the closed circles represent 
(J) and (Jin} , respectively, when D = 1. The dotted line 
and the closed squares represent (J} and (Jm} , respectively, 
when D = 5. As in Fig.Q all the quantities are dimensionless. 



and 2TR'(oj) are represented as functions of u> when 
D = 5 and a = 10. It exemplifies the violation of the 
FRR in a non-equilibrium steady state. In this case, it 
should be noted that v s = because the selected po- 
tentials are symmetric with respect to the reflection of 
the variable x. We then calculated the right-hand side 
of Eq. JSJ, which amounted to 16.891, by integrating the 
difference between C(u>) and 2TR'(u>) over the entire fre- 
quency domain. On the other hand, the rate of energy 
input for Model A has been calculated as [lil 



(•^in)o 



a I [Uiix) 



Uo(x)} [P$\x)-P[\x)] dx, 



where (P^ix), 



Pf(x)) 



is the stationary 



(8) 

solution of 



Eq. J7J). As a result of the energy balance, the rate of en- 
ergy input (Jin)o should coincide with the rate of energy 
dissipation (J) irrespective of the validity of Eq. JSJl. 
Subsequently, if ( Ji n ) coincides with the right-hand side 
of Eq. JSJ , it indicates the validity of Eq. © . Using the 
above parameters, we obtained (Jm) = 16.891. Thus, in 
this case, Eq. © has been confirmed. 

In order to demonstrate that this result is not acciden- 
tal, both the quantities defined in Eqs. JHJ and JHJ are 
represented as functions of a for two values of D in Fig. [2 
This figure clearly shows a good agreement between ( J) Q 
and (</in) ! independent of the values of the model pa- 
rameters. In this manner, Eq. JSJ has been numerically 
verified for Model A. 

Next, we provide a mathematical proof of the theo- 
rem. For simplicity, we restrict our argument to Model 
B in this proof. In the case of Model B, it has recently 
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been demonstrated that the response equation can be 
derived from the Langevin equation by transforming 
the force F(x(t)) as 



F(x(t)) = jv a 



K (t 



K ( s) + £/ P( s )] ds 



K ± (t~s;x( S ))-[as)+ef p (s)}ds, (9) 



where • denotes the Ito multiplication ^(|. We define 
K Q (t) = and K ± (t; x(s)) = for t < 0. This decompo- 
sition of force is determined by defining Kj_(t; x(s)) such 
that it satisfies the property 



(K ± (t;x(s))) = 



(10) 



for arbitrary t > and s. By taking the average of Eq. (JJJ 
after substitution of Eq. @, we obtain Eq. © by the use 
of Eq. (fTTTfl . Then, it is established that 



7-R(w) = K (w) + l. 



(11) 



See Ref. 20] for further details. 

We now set e = 0. Hereafter, in order to dis- 
tinctly treat the coexistence of two types of multipli- 
cation o and •, we discretize time as t n = nAt. Fur- 
thermore, we introduce a Wiener process 

E3, Wit), m 

place of the noise, £(t). Accordingly, the following no- 
tations are employed: x n = x(t n ), Ax n = x n+ i — x n , 
and AW n = W(t n+1 ) - W(t„). By the definition of 
the Wiener process, the relations, (AW n ) = and 
{AW n AW m ) = S nm At, hold. 

Combining Eqs. (JIJ and (J3J and taking into account 
the definition of the multiplication o, we obtain 



J{t n )At = F n Ax n + 0(At 3 / 2 ) 



(12) 



where F n = [F(x n ) + F(x n +i)]/2. On the other hand, 
integration of Eq. Q from t n to t n+ \ gives 



7 Ax„ = F n At + y^TAW n + <3(At 2 ). (13) 
A straightforward calculation then gives 



<J r (*n)>C 



' Air r , 



7^ s +7 



At 

2T (f »AW„) 
At 




+ 0(At 1 / 2 ). (14) 



In the limit of At — ► 0, the second and third terms on 
the right-hand side of Eq. (|14fl can be transformed as 




Next, since the discrete expression of Eq. © becomes 

oo 

F(x n ) = 7« s 

k=l 

oo 

+ ^jf^K ± (t k ,x n - k )AW n - k , (16) 



k=l 



it follows from Eq. ltTU|l that 



(F n AW n ) 



-K (At)At, 



(17) 



In addition, using Eq. 111|) . Fourier's integra- 
tion theorem, lim At ^o[^o(-At) + K {At)}/2 = 
K' {us)duj/2ir, leads to 



lim K (At) 



Substituting Eqs. (fTHft . (fTTjl and (JTHJ into Eq. JT3J), we 

obtain the theorem. 

We now present several discussions on the generality 
of the equality proved here. First, although this proof 
has been restricted to Model B, we can prove Eq. © in 
almost the same manner for Model A, by determining 
the decomposition of the force as expressed in Eq. 10. 
Furthermore, our result represented in Eq. JSJ can be 
generalized for a larger class of Brownian ratchet mod- 
els 18], including models with an inertia effect, a time- 
periodic potential, and spatially inhomogeneous temper- 
ature profile. We will provide a detailed description of 
these derivations in another paper. 

Second, we can also analyze systems with many degrees 
of freedom, such as colloidal dispersions, in a parallel 
way. As a simple example, we consider a system where N 
spherical colloidal particles in a three dimensional aque- 
ous solution are driven by an constant external force fe x . 
(See Ref. [2l| for a related experimental system.) Let 
us denote the coordinates of the particles by T = {xi} 
(i = 1, ■■■,3N), where r M = {x^-2, x 3li -i, x 3fJl ) repre- 
sents the position of the fi-th particle (/i = 1, •••,JV). 
Then, a widely used model describing the motion of the 
particles is provided as [13 



(19) 



where F^T) = Y%=i /<W-2 ~ E^=i U i r v) - 
dxi Hti,u=i c/ i" t (l T V~ T V I)/ 2 represents single-body forces 
and two-body interactions, and the noise satisfies that 
(&(*)£i( s )) = 2-fT6 i3 5(t - s). For this model, the en- 
ergy dissipated into the solvent is expressed as J(t)dt = 

Ei=i [7*i(*) - &(*)] ° Hi We can then p r ° ve 
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(20) 
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where Cjj(i) = ([ii(t) - (ii) o ][ij(0) - (xj) ]) Q are the 
cross correlations of velocity fluctuations, and Rijit) are 
the cross response functions defined as 



{ii{t)) e = {±i) +e 



3N 

E 

3=1 



R iJ (t~s)ff( s )d s + 0(e 2 ). 

J 

(21) 

The derivation of Eq. (|20|) will be presented in another 
paper, although it is a straightforward extension of the 
above argument for Model B. 

In addition, we present another generalization of 
Eq. (J5J. For Model B, by a similar argument with the 
above proof, we can derive an expression for the sym- 
metrized time correlation between the velocity and force 
I(t) = ([x(t) o F(x(0)) + i(Q) o F(x(t))}) /2 as 



lit) 
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C(u) - 2TR'{u>) 



2^ 



(22) 



In case of equilibrium, this equality leads to the FRR 
[Eq. ©] since we can derive I(t) = and v s = from 
the detailed balance condition. Thus, in order to further 
investigate the physical meanings of the FRR violation, 
it might be important to study the function I(t). 

Before we conclude, let us discuss the physical signif- 
icance of the theorem. From an experimental point of 
view, it has been difficult to estimate the amount of the 
rate of energy dissipation in systems under investigation. 
The virtue of the expression in Eq. JSJ is that it enables 
us to determine the energy dissipation from experimen- 
tally accessible quantities alone, without knowing the ev- 
ery detail of the system such as the profile of the force, 
F(x,t). It is expected that the present result is also use- 
ful in the study of biological molecular motors. 

In conclusion, we have presented an equality between 
the rate of energy dissipation and the extent of the FRR 
violation for a class of Langevin systems under non- 
equilibrium conditions. To the best of our knowledge, 
no previous study has addressed this equality, with the 
exception of a similar, but not precise, expression that 
was conjectured by one of the authors 24]. However, 
we should state that an inequality between the rate of 
entropy production and the FRR violation was demon- 
strated in Ref. 25]. 

In the present Letter, we restricted our argument to 
Langevin dynamics. To extend our result for more gen- 
eral non-equilibrium systems that are not described by 
Langevin dynamics is an important open problem. Re- 
lated to this problem, one might be able to prove the 
present result on the basis of microscopic dynamics, e.g., 
the Liouville equation, at least with focusing on the vicin- 
ity of equilibrium. Experimental examinations on various 
non-equilibrium systems are also of great importance. 
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